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Abstract 

Writing (R 2 N ) = AN 2u (l + BN~ Al + CN^ 1 + ...) for the mean square end- 
to-end length (R%) of a self-avoiding polymer chain of iV links, we have 
calculated Ai for the two-dimensional continuum case from a new finite per- 
turbation method based on the ground state of Edwards self consistent so- 
lution which predicts the (exact) v = 3/4 exponent. This calculation yields 
Ai = 1/2. A finite size scaling analysis of data generated for the continuum 
using a biased sampling Monte Carlo algorithm supports this value, as does 

a re-analysis of exact data for two-dimensional lattices. 
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A polymer chain is self-avoiding due to the excluded volume effect between monomer 
units which causes an expansion or 'swelling' of the chain when compared to the free random 
walk. The central quantity of interest is therefore the mean square end-to-end length (R%). 
This is believed to have the form 

(R%) = AN 2u (l + BN- Al + CN- 1 + ...), (1) 

where N is the number of chain links, v is the leading scaling exponent, A, B, C are excluded 
volume dependent coefficients and Ai is the leading correction-to-scaling exponent. It is 
now firmly established 0-fl that in two dimensions (2D) v = 3/4 is exact. Despite this, 
there is very little agreement about the value of A x . Nienhuis |l| predicts A x = 3/2, 
while Rapaport || has argued that there is no need for a correction term other than the 
analytic correction, i.e., Ai = 1. However, many numerical studies have disagreed with 
these results, with estimates for Ai of 1.2 ||, 0.84 [[F] and 0.65 [||,(|||. These numerical 
estimates are based on results obtained from self-avoiding walks on 2D lattices. With the 
exception of a very few authors [[|,[TU| (these studies however were not concerned with the 
correction to scaling terms) it appears that little work has been done in the continuum. 
Theoretical results are also in disagreement. Besides Nienhuis's prediction, which relies on 
a mapping to an exactly solvable solid-on-solid model on the honeycomb lattice, Baker et 
al [|ll|] predict A x = 1.18 using RG arguments, while Saleur [|l2j predicts Ai = 11/16 by 



conformal invariance. Interestingly, Saleur also gives evidence for a term Ai = 1/2, but he 
then rejects this result. Perturbation expansion techniques [|13|,14|, which start from the free 



random-walk solution, have also been used to predict (R%), but these methods have resulted 
in series which are divergent in iV and v, the excluded volume parameter, and hence a value 
for Ai cannot be predicted. The obvious confusion in both the numerical and theoretical 
estimates for A 1; lack of corresponding data for the continuum, and the possibility of using 
a better perturbation expansion to determine (R%) form the motivation of this study. 

We have used a new perturbation method, which unlike previous studies, starts from a 
ground state that already correctly predicts the exact large N behaviour in 2D, namely the 



Edwards self-consistent solution |15[ . Although it has been shown |16| that the Edwards 
solution cannot be the correct form for the self-avoiding random walk end-to-end distri- 
bution function ||17[ , it has mathematically convenient features that enable a perturbation 
expansion to be performed. We believe its use here underpins the essential physics and that 
Ai thus obtained may well be exact in 2D. In path integral representation [Tj^, the exact 
distribution function, or Green's function, for the end-to-end distance R is 

/■r(L)=R / 1 r L f) r ( c) v r L r L \ 

G(R,L)= D[r]exp --/ ds{-^) 2 - - ds ds'5 2 [r(s) - r(s')} , (2) 

Jr(0)=O \ I J OS l z JO Js ) 

where L is the total chain length L = Nl, I being the step length of one link, and v is the 
excluded volume. Two problems arise in dealing with this intractable path integral. Firstly, 
divergences appear in the calculation which must be handled carefully, and secondly, the 
resulting series expansion is a power series of increasing L and v. This leads to a divergent 
result unless the value of v is assumed to be very small. This divergent property is the 
hallmark of modern critical phenomena theory whose resolution was offered by the renor- 
malization group approach p9| , p0| . Historically Edwards avoided the divergence problems 
of such an approach by replacing the point contact potential by a self-consistent field W (r) 
which in 2D is equal to vp(r)/l, where p(r) is the one-particle potential proportional to 
r~ 2 / 3 fB!5|,ETtl. Therefore 



W(r) = Cv 2 l\- 2 I\ (3) 

where C = (\/3/47r/) 2 / 3 . Thus the Edwards Green's function G_b(R, L) becomes 

r W =R „ r n ( 1 r L , fdr(s)\ 2 r L 



Ge<r - l) = L- Dlr] exp ^-7 L da "/„ w{s)ds ) ■ (4) 

Our approach relies on obtaining a better first order perturbation expansion by starting 
from the Edwards ground state and then perturbing this with the difference between the 
self-consistent field and the true point contact potential ||22j| . Thus 



rr(L)=R I j rL fdr(s)\ 2 t L 

x exp (j^ W(s)ds - ^ £ ds £ ds'5 2 [r(s) - r(s')]^ , (5) 



where the difference potential in the second exponential term is now being treated as a 
perturbation. The Fourier transform of Eq. (|5|) thus becomes 

G(k, L) = G E (k, L) + G 1 (k, L) + G 2 (k, L) + (6) 

where Gi(k, L) and ^(k, L) are the first order terms in the perturbation expansion. By 



using the method of Fourier and Laplace transformation as in ||14|| , we derive the following 
functions to leading order in L 

G E (k,L) = l-^(A^ + ^L + ...\ 

G 1 (k,L) = ^k 2 A 2 L 3 / 2 + ..., (7) 

G 2 (k, L) = 2L$ - 2* - k 2 (^-^L 5 / 2 + ^-L 2 - A^/L 3/2 + ...), 
K 1 v 5 AB ' 

where A and B are excluded volume and step length dependent quantities that appear in 
the 2D Edwards solution, and $ and which also depend on A and B, are well behaved 



convergent integrals in the large L limit |HJ . The calculation of these functions required the 
exact form of the L dependent normaliztion factor of the 2D Edwards solution. However, 
unlike the normalization for the free walk solution which leads to logarithmic divergences 



and hence the introduction of a cutoff e, see [14], the normalization for the 2D Edwards 
solution has a form such that a term corresponding to e appears naturally. As we know 
its exact L dependence, integrals which would otherwise diverge remain controlled. The 
subsequent calculation of (R 2 ) from these results gives 

It should be noted that in deriving this result, unlike previous perturbation calculations, 
no divergences are encountered, and no restriction is placed on the value of v since the 
series is convergent in L. When comparing Eq. (§) with Eq. ([[]), we see that it predicts a 
value Ai = 1/2, as well as a negative value of the coefficient C. If we assume a/l = 0.5, 
corresponding to the maximum excluded volume v = ira 2 , we calculate A ~ 0.793, B ~ 1.13, 
$ w 0.107 and $ w 0.226. Substituting these into Eq. (§) the scaling amplitudes of Eq. flU) 



become A pa 0.50, B pa 2.65 and C ~ —2.07. However these 'mean field' values for the 



amplitudes are not expected to agree well with numerical or exact results [[23[ . 

We now turn to numerical studies. In order to create 2D self-avoiding chains in the 
continuum we have used a biased sampling Monte Carlo method dating back to Rosenbluth 



and Rosenbluth EM. Although more efficient algorithms exist for creating longer chains 



25j we are unsure about their reliability for studying the correction-to-scaling terms in the 



continuum. As our chains are in the continuum, the simulation procedure is considerably 



more complicated than that in ||24j| . Given a chain consisting of n circles, the area available to 
the (n + l)th circle must be determined. If no areas are large enough, the chain is discarded 
and a new one started, otherwise the position of the next circle is picked randomly from the 
available areas and the chain at that position weighted with a factor #/(27r — /?), where 9 is 
the total available angle and (3 is the angle excluded by the (n — l)th circle as no doubling 
back is permitted. This weighting factor distinguishes polymer statistics from the 'true' 
self-avoiding walk pB| . In this way relatively long chains can easily be built even with the 
maximum excluded volume. We have checked the program against the analytic solution for 
the three step walk and against a similar simple sampling program up to much longer 
chain lengths, with satisfactory results. The method has also been tested against Guttman's 
exact ennumeration data ||. 

To analyse our data, we have used the finite size scaling method of Privman and Fisher 
28] which is based on the cancellation of leading terms. We plot the estimating function 



' " iV A — (N — ]fc)A • (9) 

Assuming Eq. (P with ^ = 3/4, then 

A N>1 (A) = A + (1 - l/A^ACN- 1 + (10) 

when A = Ai. The curves for different N will thus cross at a point close to the correct 
value of A and Ai assuming \C\ is small compared to \A\. The insert of Fig. [I] shows 
this technique applied to the exact square lattice data of Guttmann || for values N = 15 



to N = 27. In this case k = 2 is used and the resulting Ajv i2 (A) data averaged with 
previous v4jv_ 1j2 (A) result to eliminate the odd-even effect. As can be seen the curves cross 
at a value A « 0.65 and A « 0.765, in agreement with Ishinabe [[J, while Ai = 0.65 
was also reported by Privman M using the same technique on the triangular lattice data 



of Grassberger However, these authors assumed \C\ -C \A\, even though no clear 

evidence was given to support this assumption. Figure [l] also shows the results obtained 
using the same method on data calculated from Eq. ([!]) (hereafter called 'simulated data') 
with A = 0.760, B = 0.227, C = 0.18 and Ai = 0.5. These values of A, B and C were 
obtained from a least squares fit to the data || having set Ai = 0.5. The resulting curves 
are virtually indistinguishable from the exact data even though the values of A x differ. We 
have confirmed that any positive value of C will shift the crossing point to higher values of 
A. It thus appears that for the lattice Aj could be as large as 0.66, but Aj = 0.5 is possible. 
The same procedure using exact triangular lattice data |29| yielded excellent agreement with 
simulated data for A = 0.704, B = 0.175, C = 0.128 and A l = 0.5. 

Figure |2] shows the resulting curves when the estimating function Eq. (|9|) is applied to 
our continuum data for chain lengths of N = 10 to N = 25 with the maximum excluded 
volume ratio of 0.5. The estimate of (R%) for N = 25 comes from averaging approximately 
1.87 x 10 8 walks, resulting in an error of less than 0.1%. Clearly the data show no sign 
of crossing, and we suggest that this is due to a large negative C. An estimator, similar 
to Eq. (§), but which gives simultaneous estimates of B and A x || was also studied, but 
it too showed no evidence of crossing. Due to the large C value and small random errors 
in our Monte Carlo continuum data it is difficult to use other graphical techniques 0,0] to 
determine Ai. We therefore use the above method of comparison between simulated and 
Monte Carlo data. 

We first assumed a value of A x = 0.666, in agreement with ||], and using a least squares 
curve fit to our data from N = 10 to iV = 25, we found a best fit with coefficient values 
of A = 1.023, B = 0.491 and C = —0.681. We then used these values to create simulated 
(R%) data to which we applied the Ajy^(A) vs. A analysis. The resulting curves are shown 
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in Fig. [3[ They are very different from those of the Monte Carlo continuum data. We then 
assumed a value of Ai = 0.5 from which we obtained a best fit with coefficient values of 
A = 0.990, B = 0.489 and C = —0.839. The resulting curves, as shown in Fig. |], are in 
excellent agreement with the continuum data. The large negative C should be noted. Thus 
to assume C is negligible, as has often been done when analysing 2D lattice data is 
a mistake and could lead to poor estimates of Ai. 

In conclusion, we have presented strong evidence that Ai = 1/2 for 2D chains in the 
continuum. Although this value does not agree with any of those suggested by other authors 
for 2D chains on the lattice |||],[7|J|, their data are compatible with Ai = 0.5 when the effects 
due to the next order term in Eq. ([J) are considered. Unless there is a breakdown of the 



universality of both v and Ai p0| , p0| , we suggest that Saleur's |L2| rejection of Ai = 1/2 
should be re-examined. Since our new perturbation method and our Monte Carlo analysis 
agree, we suggest that the leading correction-to-scaling term in two-dimensions is Ai = 1/2. 
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FIGURES 

FIG. 1. Plot of A Nt2 (A) vs. A for simulated square lattice data with A = 0.760, B = 0.227, 
C = 0.18 and Ai = 0.5. The insert shows the same plot for the exact square lattice data of Q. 

FIG. 2. Plot of A/v,i(A) vs. A for our Monte Carlo 2D continuum data. 

FIG. 3. Plot of A NA (A) vs. A for simulated continuum data with A = 1.023, B = 0.491, 
C = -0.681 and Ai = 0.666. 

FIG. 4. Plot of Ajv,i(A) vs. A for simulated continuum data with A = 0.990, B = 0.489, 
C = -0.839 and Ai = 0.5. 
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